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The generalization of the dimer model on a two-leg ladder is defined and investigated both, ana- 
lytically and numerically. For the closed system we calculate the Landauer resistance analytically 
and found the presence of the point of delocalization at the band center which is confirmed by the 
numerical calculations of the Lyapunov exponent. We calculate also analytically the localization 
length index and present the numerical investigations of the density of states (DOS). For the open 
counterpart of this model the distribution of the Wigner delay times is calculated numerically. It is 
shown how the localization-delocalization transition manifest itself in the behavior of the distribu- 
tion. 
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I. INTRODUCTION 

One of the widely accepted results of the disordered 
systems is that in one-dimensional spaces all the elec- 
tronic states in independent site-energy random (diago- 
nal disorder) Anderson model with non-random hopping 
amplitudes are exponentially localized The diffusion 
or long-range transport of the initially localized particle 
in such systems is absent regardless the strength of the 
disorder. The same result is proved rigorously for quasi- 
one-dimensional systems 0. 

On the other hand, it is well established that the 
type of disorder can have a strong effect on localiza- 
tion properties of disordered systems. For instance, in 
one-dimensional systems with off-diagonal disorder (ran- 
dom hoppingmodels) delocalization takes place at the 
band center In quasi-one-dimensional systems with 
off-diagonal disorder the picture is more complicated. It 
appeared 0, Q , that in a weakly disordered N-\eg quasi- 
one-dimensional disordered tight-binding hopping model 
there is a delocalization transition at the critical energy 
Ecrit — if and only if N is odd (including N — 1 case). 
The wave-functions remain localized for even N . The 
localization properties of two- and three-channel tight- 
binding Anderson model with nearest-neighbor inter- 
chain hopping was studied in Q . 

Another possibility to have localization-delocalization 
transition in one dimension is the correlated disorder, 
when the random variables are correlated at short or long 
distances. There is intensive current interest both the- 
oretical (seeHJ| IIEIIIEHIiIEEEH 

f l20l I21I I22I |23j and Refs. therein) and experimental 
25] concerning analyze of conditions, under which 
delocalized states can appear. This kind of models with 
correlated disorder give explanation of high conductivity 
of polymers such as doped poly-aniline and also describe 
the transport properties of random semiconductor super- 
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lattices 0, 0, 0] • Some transport properties of many- 
mode waveguides with rough surfaces were studied in [2|| 
and the possibility of perfect transmission due to some 
specific long-range correlations in the surface profiles was 
investigated. The scaling properties of the models with 
correlated disorder were studied in [23, l2ll |23| . In [2(j 
it was shown, that universality of one parameter scaling 
breaks down in presence of long range correlated disor- 
der, for instance in the periodic-on-average systems. It 
is necessary to point out, that the correlated disordered 
systems do not belong to Zirnbauer's [2?J and Caselle's 
[23 classification of fully disordered systems on the basis 
of Cartan's classification of symmetric spaces. That hap- 
pens because, due to correlations, the Brownian motion 
of the effective degrees of freedom carries out on the re- 
stricted part of the coset space. This restriction is general 
characteristic of all correlated disordered systems. 

So far the discussion of correlated disorder was applied 
only to strictly one-dimensional systems. In order to ana- 
lyze the existence of localization-delocalization transition 
in higher dimensional lattices with correlated disorder, it 
is natural to make the first step in this direction — to 
study the quasi-one-dimensional systems. To investigate 
the existence of critical point and address the problem 
of delocalization transition in this case analytically is 
an important task. The exact calculation of Landauer 
localization length and critical indices will allow us to 
understand whether numerical simulation of Lyapunov 
exponent shows truly critical point (where the localiza- 
tion length is infinite) or an extended state (where the 
localization length is finite, but exceeds the system size). 

In this paper we will consider the extension of the ran- 
dom dimer model (RDM) 0, S, S E3, E3, E2, O S 
to quasi-one-dimensional case and show the existence of 
delocalization transition for particular realization of cor- 
relation. Namely, we consider a two-leg ladder chain (see 
Fig.l) and a tight-binding model of random binary al- 
loy with on-site potential fields taking values e a , €b from 
the group I2 , randomly assigned to the lattice sites with 
probabilities p and (1— p), correspondingly. Here we will 
focus on a particular case, when the same site poten- 
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FIG. 1: Two leg ladder 



tials appear always by quartets disposed on four nearest 
neighbor sites creating squares, as it is marked on the 
Fig.l. This is most natural generalization of the RDM in 
case of ladder models, and one can call it random ladder 
dimer model (RLDM). 

The Hamiltonian of the RLDM is given by 

N 

rj=l;i=l,2 

^f C n,'i C »V'+l ~1" e n,i c n ,i C n,i} ■ (1) 

Here t v and th are constant vertical and horizontal hop- 
ping parameters correspondingly, e n j is the random ex- 
ternal potential at site n of the i-th chain. N de- 
notes the number of atoms in each chain. We impose 
a constraint on the randomness by demanding that four 
£n,i (n = 2k — 1,2k; i = 1,2) around the plaquette (see 
Fig.l) always have the same, but randomly chosen value 
from the set {±m} 



e2k-i,i — £2fc-l,2 — £2fe,l — £2fe,2 € {±m} , 

N 
~2' 



1. 



(2) 



In section |n] we study the closed system described by 
Hamiltonian fl}. We use the technique introduced in 
PI, E3> and further applied to the RDM in 
in order to investigate the RLDM. We calculate analyti- 
cally the dimensionless Landauer exponent of the model 
and compare it with numerical simulations of the Lya- 
punov exponent (Fig. 2 and Fig. 3), calculated iteratively 
in standard way |29j, |30j. The exact calculation of the 
Landauer resistance shows the existence of one real crit- 
ical point with critical energy E cr n = in the case of 
m = t v /th- Otherwise all states are localized. 

In the case of m = t v /th < 1 we analytically analyze 
the divergence of the inverse of Landauer exponent and 
find, that there is delocalization with the critical index 
v = 2. On Fig. 2 and Fig. 3 we see very good coincidence 
of the half of Landauer exponent with the Lyapunov one 
around the critical point. 

The critical behavior is essentially different when 
m = t v /th = 1. In this case also the localization length 
diverges at the point E cr u = 0, but with index v = 1/2. 
The comparison of Landauer and Lyapunov exponents 
for this case is presented on Fig. 4. 

The numerical analyze of the density of states (DOS) 
shows anomalous behavior at the band center in the case, 



when t v = th (see Fig. 5) and regular behavior for t v < th 
(see Fig. 6). 

Section ITTll is devoted to the open counterpart of the 
system described by Q. There we show that far from the 
delocalization point the distributions of the Wigner de- 
lay times are very similar to those found for the random 
models with uncorrelated disorder (Fig. [Hand Fig. IHJ • At 
the critical point the behavior of the distribution depends 
on m. Namely, for m = 1 the distribution of the re-scaled 
Wigner delay times tends to the S- function (Fig. 0) , cor- 
responding to the deterministic, ballistic propagation. 
While for m < 1 this distribution is not deterministic, 
but still bounded (Fig. [113} , indicating the ballistic prop- 
agation with velocities distributed in some interval. 



II. LANDAUER RESISTANCE, CRITICAL 
EXPONENTS AND DENSITY OF STATES 

For one particle eigenstates ip n (E) of the eigenenergy 
E the Schrodinger equation of the Hamiltonian 0} has 
the following recursion form: 

th{lpn+l,i + tpn-l,i) + t v tpn,i+l = {E — t n ,i)^n,i- (3) 

Now we can rewrite (01 in the matrix form 




(4) 



by introducing 4x4 transfer matrix T n as 
E 

T n = 



-e„,i ~t -1 

-t E-e n . 2 -1 

1 

10 



(5) 



Here we re-scaled the energies and the vertical hopping 
parameter t v (introducing t instead) by the horizontal 
hopping parameter th- One can easily find out, that the 
4x4 transfer matrix T n possesses the condition T+ JT n — 
J with 



/0 -i 

-i 

1 
\0 i 



(6) 



which means that T n is an element of the group SU (2, 2). 
If we introduce a new matrix M n = Y\ i=n Ti, then we can 
relate the wave functions (4> n +i,i, VVi+1,2, ipn,i> VV2) with 
ones at sites and 1 of both chains 





(7) 



and Mn will denote the total transfer matrix of the sys- 
tem. From the definition of Mjv it is obvious, that it 
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has the same property M^JMjv = J- This condition 
is equivalent to the flux conservation in the theory (i.e., 
Hermiticity of H). 

In a physics of disorder particular interests has the 
limit of eigenvalues of the matrix ln(MjvMy) w when 
N — ► oo. According to Oseledec theorem [^lj this limit 
exists and the maximal eigenvalue is called Lyapunov ex- 
ponent. Due to self averaging property it is determined 
by the maximal eigenvalue of the matrix M/v 



y[e] 



iLy 



N 



lim (In || M 



N 



(8) 



where ||M W || 2 = EaA MN ^( M N^ = Tr(M N M+) and 
the average is taken over random distribution of the ex- 
ternal potential taking into account the correlations (if 
they are present in the model under consideration). 

It is hard to calculate the Lyapunov exponent exactly 
and establish the presence of derealization point. An- 
other variable, which shares the essential characteristics 
of the localization-delocalization transition but can be 
calculated exactly is the dimensionless Landauer resis- 
tance, which is defined as the ratio of squares of the 
modules of the reflection and transmission amplitudes. 
According to arguments of Anderson et al. [32j in one- 
dimensional systems the Landauer exponent (defined in 
Eq. 10) is simply the double of Lyapunov one. The pro- 
portionality of these two exponents was confirmed around 
of the critical point in the paper [33j . Considering in most 
general terms one dimensional models with transfer ma- 
trices belonging to SU(1, 1) authors of [3^ have shown, 
that in the lowest order of perturbation expansion with 
some parameter measuring the distance from the criti- 
cal point, the Landauer exponent is twice the Lyapunov 
one, just as predicted by |32j| . This means, that critical 
indices of both exponents are equal. In quasi-one dimen- 
sional models, which have transfer matrix from SU (n, n), 
again one can expect the coincidence of critical indices. 

The exponential increase of Landauer resistance with 
distance is determined by the Landauer exponent , which 
can be defined as 



j(E) = lim ln(|| MjvHK 



jv 



(9) 



In the articles ^ was shown, that the Lan- 

dauer resistance (and corresponding exponent) can be 
calculated exactly just by reducing the direct product 
Mjv <8> Mjy of fundamental representations of transfer 
matrices M of SU(1, 1) group to the adjoint one. We 
apply this technique here for the group SU{2 1 2). In or- 
der to calculate this direct product exactly, we use the 
known representation of the permutation operator via 
generators t m {h — 1, . . . , 15) of the sl(4) algebra as 
P = j (I + r p <g) r* 4 ). In the matrix elements this for- 
mula looks like 

%lS%=\l^+(T»)% 1 (T lt )%], (10) 

where we suppose summation by the repeating indices 
ji. Among of generators r M there is one, which coincides 




FIG. 2: The Lyapunov (curve B) and the half of Landauer 
(curve A) exponents versus energy in the whole interval 
[-2.5,2.5] of existing eigenstates in case of m = t v /th=l/2. 
The chains of the length N=40000 at 60 values of energies 
were considered. 




FIG. 3: The Lyapunov (curve B) and the half of Landauer 
(curve A) exponents in the vicinity of the critical point E=0 
versus energy E in case of m = t v /th = 1/2. The chains of 
the length N =40000 for 60 values of energies in the interval 
[-1.1, 1.1] were considered 



with the metric J defined in 0. We denote the corre- 
sponding index [i as J, namely t j = J. 

By simple multiplication of (|10fl by Tj and from 
the left and right hand sides correspondingly one can 
express the direct product of Tj and Tf via their adjoint 
representation: 



(t 3 TAt+% = ~{J) a {JZ + ~(^J)£,Af (Jt^.(ii) 

Here the adjoint representation A„ of T n is defined by 



0' 



Ar = -Tv{T nT »T+T») 



(12) 



and is an 15 x 15 matrix depending on the parameters of 
the model at the site n of both chains. 

Now, in order to calculate the average of Landauer ex- 
ponent, we use l|ll|) and decompose the direct product 
of two transfer matrices corresponding to the whole sys- 
tem of 2A sites. In the model under consideration we 
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attach the same site potential e to the quartet of sites 
(see Fig.l) and calculate the average over random distri- 
bution of values e a and ej, with probabilities p and 1 — p 
respectively. Therefore we should consider the A 2 as a 
constituent transfer matrix and average it in a product 



1 1 ( N/2 \ 

(M N M+) = —J(g)J+ — (r M J) ® {Jt v ) I II< A i>J 

It is clear, that 

(A 2 )^ = \ Tr (TjV[T/]- V) 



and the average is 



A = (A 2 }=M 2 „M + (1-p)A 2 M, 



//// 

(13) 
(14) 

(15) 



where Aj(e ) and Aj(e;,) are calculated for the site po- 
tentials e a and tb respectively. The matrix elements of 
A can be found easily from Ijl4|l and Ijl5jl . Hereafter we 
take e a = —to, ei, = +m without lose of generality and 
p = 1/2 for simplicity. 

From the formula l|13|l one can easily obtain 



M || 2 ) =4(A W / 2 ) JJ 



x N/2 
max ' 



(16) 



where the biggest eigenvalue X max defines the asymp- 
totic behavior of Landauer resistance with the distance. 
In (|16|) superscript J J means J J element of the matrix 
k N f 2 (see notion after formula (|10p *> . 

According to © the Landauer exponent j(E) con- 
nected with the maximal eigenvalue as 



j(E) = - log[A max ] 



(17) 



and defines the inverse of Landauer localization length 
£(E) = lh(E). 

The delocalized states correspond to the critical points 
where £(E C ) — oo, or X max — 1, and can be found by the 
condition det \ I — A |= 0. In our case this characteristic 
equation is given as 

m 4 l{E-t) 2 -m 2 ][(E+t) 2 -m 2 ][m 4 -E 2 (-l+m 2 )}t 2 = 0. 

From this equation one can find following candidates for 
critical points 



Eic 
E-ic 
E-^c 



t ± TO, 

—t ± m, 
TO 4 /(m 2 - 1). 



(19) 
(20) 
(21) 



But since the Landauer exponent defined by the max- 
imum eigenvalue of A one need to find out whether this 
points belong to maximum eigenvalue or not. The full 
characteristic equation det | AI — A |= for eigenval- 
ues A factorizes into the product of three polynomials 
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FIG. 4: The Lyapunov (curve B) and Landauer (curve A) 
exponents in close vicinity of the critical point E=0 versus 
energy E in case of m = t v /th = 1. The chains of the length 
N=40000 for 60 values of energies in the interval [-0.1, 0.1] 
were considered 

p[E] 




FIG. 5: The density of states dependence on energy for the 
boundary case m = e a — eb = 2 calculated numerically for the 
ensemble of 4000 coupled chains with N = 50 sites in each 
chain. 



Ui(X, E, t, m)J7a(A, E, t, m) (Ua(X, E, t, to)) and gives us 

Ui(X,E,t,m) = (1 - 2to 2 ) 2 + [(2to 2 - 1)(3- 
2to 2 + (2E + E 2 - to 2 - 2(1 + E)t + t 2 )(E 2 - 
to 2 - 2E(1 +t) + t(2 + t)))] A + [3 + E 4 + 
to 4 - AEH - At 2 + t 4 - 4M(3to 2 + t 2 - 2)+ 



m 



2 {6t 2 - 4) + E 2 (6m 2 + 6t 2 - 4)] A 2 — A 3 = 



(22) 



U 2 {X,E,t,m) = (1 - 2to 2 ) 2 + [(2to 2 - 1)(3- 
2to 2 + (2E + E 2 - to 2 + 2(1 + E)t + t 2 )(E 2 - 
m 2 - 2E{1 -t)- t(2 - t)))] A + [3 + E 4 + 
to 4 + 4EH - 4i 2 + t 4 + AEt{im 2 + t 2 - 2)+ 
m 2 (6t 2 - 4) + £; 2 (6to 2 + 6< 2 - 4)] A 2 - A 3 = 0, 



(23) 



LT 3 (A, E, t,m) = -(l- 2m 2 ) 2 - (2to 2 - 1) [E 4 + m 4 + 
(t 2 - 2) 2 - 2E 2 (2 + to 2 + t 2 ) + m 2 (6t 2 - 4)] A + 
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FIG. 6: The density of states dependence on energy for the 
case m — (e a — £b)/2 = 0.5 calculated numerically for the 
ensemble of 2000 coupled chains with N = 50 sites in each 
chain. 



2 [6?7 
3m 2 



3 - 3m 4 



E 4 (m 2 



2(m 4 



- 2)t 2 + (m 2 - l)t 4 - 2E 2 (m 4 - 2 + M 2 - 
3m 2 (t 2 - 1))] A 2 + [E 4 + m 4 + (t 2 - 2) 2 - 

2)] A 3 — A 4 = 0. 

(24) 



2m 2 (2 + t 2 ) + 2£ 2 (3m 2 - 1 2 



Easy to see, that U2(\,E,t,m) = Ui(X,E,—t,m) 
The critical points E\ Cl E2c and E^ c are zeros of 
U\(\, E, t, m), U2(X,E,t,m) and ^(A, E, t, m) respec- 
tively. 

Cubic and quartic equations 1(22124(1 can be solved an- 
alytically. It appears, that the roots of quartic equa- 
tion 124(1 are either less of one or complex valued, there- 
fore they are out of interest. Roots of cubic equation 
f/2(A, E, t, m) = can be obtained from <7i(A, E, t, m) — 
by simple change t — > —t. Both of them behave as 
roots of the corresponding equation in one-channel RDM 
(see formula (14) and Fig.l in [l^ ) but are shifted on ±t 
respectively. Therefore, in general situation (t ^ m), the 
roots of one characteristic polynomial arc in a shadow of 
the brunch of the other polynomial and the maximum 
eigenvalue is always more than one. But if one will fine 
tune m = t , then from the expressions 1(19120(1 of E\ c 
and i?2c it follows, that the point E c — is common 
critical point and, therefore, it is a critical point of max- 
imum eigenvalue X ma x- In the limiting case t — ► 0, when 
the transverse hopping is absent and we have two decou- 
pled dimer chains, there is no any shift of two curves, 
they just coincide. Then we obtain the same dereal- 
ization points ±m, defined by formulas (|19I2U|) . as in 
one-channel RDM. We find also, as in one-channel RDM 
0, H S 13 113 j the presence of critical point only for 
case of 77i < 1, when each of branches can have critical 
points (H3 and (|2T)|) . 

The presence of derealization point E = can be 
proved rigorously analyzing the Lyapunov exponent. We 
have two transfer matrices (depending on ±m) randomly 
distributed in a product along a chain. The maximal 



eigenvalue of this product after statistical averaging just 
defines the Lyapunov exponent JSJl of the model. If there 
is a particular value of E, where this two transfer ma- 
trices commute (this means all matrices in the prod- 
uct commute) then the eigenvalues of the product are 
given by the product of the corresponding eigenvalues of 
the constituent random matrices. If, at the same time, 
the eigenvalues of both transfer matrices at this point 
are laying on the unit circle (modules of eigenvalues are 
equal to one), we have a point of derealization. In or- 
der to see that, let us diagonalize the transfer matrix 
Mjv by the matrix S, which is also diagonalizing ma- 
trix of the individual Ti 's along a chain. One obtains 
iLyapunov = Hm^oo m \nTr(S MdS' 1 (5" x )+ M+ S+) , 
where the eigenvalues in Md, as a product of individual 
eigenvalues of X^, also have module equal to one. There- 
fore, the argument of logarithm is bounded from above 
by some constant independent of N and the limit N — > oo 
is zero. It remains now to check, that indeed the squares 
of transfer matrices © for t n> i — ±m are commuting at 
E = and for t = m. For t < 1 they have two eigenval- 
ues equal to —1 and two complex conjugate eigenvalues 
on the unit circle (e ±ta ). This prove the presence of the 
derealization point. 

In Fig. 2. we present the graphic of the half of Landauer 
exponent j(E) versus E (curve A), defined by the max- 
imum solution of cubic equations 122I24|1 for the generic 
value m = t v jth — 1/2- We present here also the Lya- 
punov exponent 7 (curve B) obtained by the numeri- 
cal simulations of the chains of the length TV = 40000 
at 60 values of energies in the whole energy interval 
[—2.5, 2.5]. The same exponents in the vicinity of critical 
point [E 6 [—1-1, 1-1]) are presented on the Fig. 3 in or- 
der to demonstrate the coincidence of critical behaviors, 
as it expected after the articles (3^ and |3^| . The critical 
index v is 2 in this case. 

In Fig. 4 we present the Landauer and the Lyapunov 
exponents versus E for the boundary case m = t v /th = 1. 
We present results for close vicinity of E = and one 
can see exactly the same critical behavior of the both 
exponents. It is necessary to mention, that in contrast 
to the generic case m =/= 1 here the multiplication factor 
1/2 is absent. 

For m > 1 the derealization point E c = disappears. 

The critical index v of the Landauer exponent one can 
obtain by expanding the solutions of the equations l|22t - 
I23|l near the critical point E c = 



±E 2 + k 



E 3 + 0[E] 4 , if m = t = \ 



1(E) 



E - \E + 0[Ef/ 2 , if m = t = l. 



(25) 



As we see the critical indices are different for generic 
777, = 1/2, v = 2 and boundary m = 1, v — 1/2 cases. 
Fig. 3 and Fig. 4 demonstrate that difference. The sim- 
ilar situation was in the one-channel RDM [ijj. The 
index v = 1/2 is equal to the corresponding index in 
one-channel RDM in the same boundary case m = 1, 
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FIG. 7: Distribution of the Wigner delay times r for E — 0.3, 
t = 1, m = 0.1 and N = 10 4 (circles). The line represents the 
best fit by formula l|88|l with the fit parameter a. 

when one is approaching to the critical poin ts E c = ±1 
from the outside of the region [-1,1] (see[l9| ). After the 
discussion presented above this fact is natural to expect. 

In order to fill up the exact results from above, we 
have analyzed the density of states in the model by the 
numerical diagonalization of the Hamiltonian (0. 

In Fig. 5 the DOS versus energy for the boundary case 
in = e a — £fc = 2 is presented. Here we calculated the 
DOS for a system with TV = 50 sites in each chain. The 
result comprise the average over 4000 realizations of the 
disorder. Here one can clearly see the anomaly of DOS 
at band center E = (in accordance with |34(), corre- 
sponding to derealization transition. A similar anomaly 
is present in RDM at the both critical energies e a , £fc in 
the same case when c a — 65 — 2. 

In the Fig. 6 we present the normalized DOS for the 
generic case m — 0.5, calculated numerically for the en- 
semble of 2000 coupled chains of length N — 50 each. 
This case corresponds to a different type of critical be- 
havior with respect to the previous one and the DOS is 
finite in whole energy interval. 

III. WIGNER DELAY TIME 

In this section we consider the properties of the open 
counterpart of the model studied in section [H] . Namely 
we investigate the statistical distributions of the Wigner 
delay times for various regimes. The Wigner delay time is 
defined as the derivative of the total phase of the scatter- 
ing matrix with respect to the energy r = — i 9ln ■ 
This quantity was intensively studied in recent years: for 
chaotic HH, localized [M El UlI , diffusive paH3 sys- 
tems as well as for systems at criticality [SfiL l40j|. 

In order to study an open system we attach a perfect 
lead to the semi-infinite disordered sample. Thus e n ,i = 



0, t = 1 for n < and they take the same values as in 
the previous section for n > 0. Let us first assume more 
general situation of M channels and derive the expression 
for the Wigner delay time. 

The wavefunction in the lead is a superposition of the 
plane waves in the longitudinal direction and the stand- 
ing waves in the transverse direction: 

j=l { V smK J y/SW-Kj J 

(26) 

where (j>{ = sin ^ m+\ ^) 1S ^ ne transverse eigenf unc- 
tion and kj is the wave-vectors of jth mode determined 
by the equation E — 2th cos kj + 2t v cos ^ A 7+i ) ■ ^ n orc ' er 
to decouple different modes it is useful to introduce new 
variables: 

Wn = U$ n , (27) 

where W n = (W n ,i, ■ ■ ■ , W„,m) T , i/>n = 
(ip n ,i, ■ ■ ■ , ipn,M) T and U is M x M matrix with 
elements Uij — (/)*. In new variables each component 
of the wavefunction in the lead contains contribution of 
one particular mode only: 

W„, ; = - Ape'*'" + _J?L=e-*" (28) 
V sin ki v sin k\ 

In order to rewrite (QJ) in terms of the new variables we 
notice that the generalization of the transfer matrix (0 
for M channel case is 2M x 2M block-matrix: 

r« = (^o 1 )- ( 29 > 

where C n is a three-diagonal matrix generalizing the left 
upper block of the transfer matrix (J5J. Now the analog 
of (J2J in new variables takes the similar form: 

(l +, H(tJ' <-> 

where the new transfer matrix T n is giving by 

f n = ( C » - 1 V C n = UC n V- x (31) 

Let us now consider the finite sample of length N and 
impose Dirichlet boundary condition at n = iV + 1. Then 
the analog of equation 10 reads 

(6 \_(P$ Pk 2 \(7^[^+^) \ 32) 

[W N ) - { P% P$ ) [ ^= ( e -^l N + ) {? ) 

where Pjv is the total transfer matrix Pjy — Yl n =i ^ n ano - 
matrix K is diagonal K = diag(fci, . . . , fcjw)- 
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FIG. 8: Distribution of the Wigner delay times r for E — 0.3, 
m — t — 1, and N = 10 4 (circles). The dashed line represents 
1/t 2 behavior. 



FIG. 9: Distribution of the rescaled Wigner delay times f = 
t/N for E — 0, m — t — 1, and various system sizes: (a) 
N = 100, (b) N = 1000, (c) N = 10000. 



The scattering matrix for our geometry consists of only 
one non-trivial block — the reflection matrix Rn, which 
relates the incoming wave amplitudes An to the outgoing 
wave amplitudes B^: 



Bn — RjqA 



JY 



(33) 



Using this definition and relation (|32|) one can easily cal- 
culate the reflection matrix: 



Rn 



I + V 'sin KF N [VsinK 



V sin KF N (VsinJf) 



iK 



iK 



N 



pl2 

•TV 



(34) 



Taking the derivative with respect to energy from this 
expression one can find the formula for the Wigner delay 
time: 



tn 



2ImTr 



F 



N e 



-iKl 



The matrix Fn and its derivative with respect to energy 
can be calculated recursively. To this end we increase the 
length of the sample by one and see how the total transfer 
matrix changes. Then from the definition of Pjy an d 
Eq. I|31|) one can derive the following recursion relations: 



Fn+i 
9Fn + i 
dE 



-{Cn 
Fn+i 



, dF N 



dE 



F NA 



(36) 



Starting from the initial values F\ = — (Ci) _1 , 
dFi/dE = Fi and iterating relations (l.'iOII one can calcu- 
late Fn and dF^/dE for any system size N. Inserting 



these values into Eq. (|35|l we obtain the Wigner delay 
time tn- 

Below we present the results of numerical calculations 
of the distribution of tn for the random ladder dimcr 
model. In these particular case the matrix C n becomes 
diagonal: 



E-t- e n 

E + 1 - e n 



(37) 



therefore the recursion relations l|36|l are reduced to two 
decoupled equations. That makes the numerical simula- 
tions especially efficient: in all calculations we used 10 5 
realizations of the random potentials and the system size 
up to N = 10 4 . 

Let us now using the knowledge gained in the previ- 
ous section consider different regimes. We start with the 
localized regime. In Fig[7]the distribution of the Wigner 
delay times is presented for t = 1, m = 0.1 and E = 0.3. 
The distribution is described very nice by the formula: 



V(r) 



-a/r 



(38) 



which is derived for the distribution of r for a one- 
dimensional Anderson model in the weak disorder limit 
[36l l37| . This indicates that correlations in disordered 
potential though renormalize the localization length, but 
does not change the form of the distribution in the weak 
disorder limit. The increasing the strength of the disor- 
der modifies the short time behavior of "P(r) (Fig. [SJ . 
However the typical 1 /r 2 tail of the distribution remains 
unchanged. This result is again very similar to one ob- 
tained for a one-dimensional Anderson model |36| . 

Next we concentrate on the critical regime. The most 
surprising behavior of V(t) is found for t = m = 1 and 
E = 0. As it follows from Fig. the distribution of 



P(T) 
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FIG. 10: Distribution of the rescaled Wigner delay times 
f = t/N for E = 0, t = m = 0.7 and various system sizes 
N = 100, 1000, 10000. 



the rescaled Wigner delay times f = t/N goes to the 
S- function in the thermodynamic limit. This means that 
with the probability one the particle propagates through 
the sample ballistically with a fixed velocity. The bal- 
listic propagation is compatible with the derealization 
taking place at this point. The deterministic character 
of the distribution is however not trivial and can not be 
deduced from the derealization only. This becomes very 
clear when we choose t — m < 1. The distribution of 
the rescaled times (Fig. HOfl is bounded (absence of long 
tails), but not deterministic anymore. In contrast it has a 
complicated oscillating structure in the middle and high 
peaks at the edges. The difference in the behavior of the 
distributions of the Wigner delay times for t = m = 1 



and t — m < 1 is another manifestation of the differ- 
ence between the critical indices and DOS found in the 
previous section. 



IV. CONCLUSION 

We have formulated and investigated the generaliza- 
tion of random dimer model to the case of ladder chain. 
As in one-channel case we found derealization point at 
E c — with two type of critical behavior when e a — e& < 2 
and e a — eb = 2. The comparison of Landauer and Lya- 
punov exponents shows that they have the same critical 
behavior in full accordance with the articles [33 and [33| . 

We studied also the distribution of the Wigner delay 
times for this model. Our results shows that its behavior 
ranges from 1 /r 2 decay typical for the localized system to 
the deterministic one at the critical point for e a — e& = 2. 

Our analytical and numerical results can be relevant 
for photoluminescence experiments on electronic proper- 
ties of GaAs-AlGaAs superlattices [24j and for experi- 
ment on a microwave waveguide with inserted correlated 
scatterers |25| . 

At the end we would like to make a following remark. 
The dimer model with the correlated disorder defined 
here does not belong to the classification of fully disor- 
dered systems on the basis of symmetric spaces j22J |28| . 
This is because the geometry of integration space over 
random variables, due to the restriction imposed by cor- 
relations, does not coincide with geometry of symmetric 
spaces. The classification of universality classes of the 
correlated disordered systems is an open and interesting 
problem. 
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